-I— > 



On the Wiener-Hopf factorization for Levy 
processes with bounded positive jumps 



A. Kuznetsov * , X. Peng t 

Current version: August 16, 2011 



o 

bD 

^ Abstract 

We study the Wiener-Hopf factorization for Levy processes with bounded positive jumps and 
arbitrary negative jumps. Using the results from the theory of entire functions of Cartwright class 
we prove that the positive Wiener-Hopf factor can be expressed as an infinite product in terms of 
the solutions to the equation ip(z) = q, where ip is the Laplace exponent of the process. Under 
some additional regularity assumptions on the Levy measure we obtain an asymptotic expression 
for these solutions, which is important for numerical computations. In the case when the process is 
spectrally negative with bounded jumps, we derive a series representation for the scale function in 
terms of the solutions to the equation ip(z) = q. To illustrate possible applications we discuss the 
implementation of numerical algorithms and present the results of several numerical experiments. 
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O 1 Introduction 

•• Assume that we want to study the way in which one-dimensional Levy process X exits a half-line or 
• i-H a finite interval. For example, we might be interested in the first passage time across a barrier, the 
^ overshoot /undershoot at the first passage time, the last time that the process was closest to the barrier, 
03 the location of the process at this time, etc. These questions are usually referred to as "exit problems" in 
the literature, and they have stimulated a lot of research in recent years due to numerous applications in 
such diverse areas as actuarial mathematics, mathematical finance, queueing theory and optimal control. 

Let us denote the supremum S t = sup{X s : < s < t} and infimum I t = inf{A s : < s < t}, and 
let e(q) be an exponentially distributed random variable with parameter q > 0, which is independent 
of the process X. It is an established fact that exit problems are closely related to the Wiener-Hopf 
factorization, which studies the distribution of S e t q -\ and I e (q)- For example, if we know the positive 
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Wiener- Hopf factor (which is defined as the Laplace transform of S e ( q )), then through the Pecherskii- 
Rogozin identity [31] we know the joint Laplace transform of the first passage time and the overshoot. 
The bad news is that for general Levy processes the Wiener-Hopf factors cannot be obtained in closed 
form, therefore the best that we can do is to try to find rich enough families of Levy processes with 
special analytical properties, for which we can say something useful about the distribution of 5 e ( ? ) and 

Let us look at the existing examples of Levy processes for which one can identify the Wiener-Hopf 
factors and the distribution of extrema. In the general case, when the process has jumps of both 
sides, this list includes processes with jumps having rational transform [29, 32] and recently introduced 
meromorphic processes [21]. The first class includes processes with hyper-exponential [7, 14, 17] and 
phase-type jumps [1, 2], while meromorphic processes include Lamperti-stable processes [4, 5, 8, 30], 
hypergeometric processes [6, 22, 26], /3-processes [19] and ^-processes [20]. In the simpler case when the 
process is spectrally negative (which means essentially that it has only negative jumps) it turns out that 
both of the same two classes provide analytically tractable formulas, however in this case there also exist 
other interesting families, such as the processes constructed in [13] (see also [23]). 

One may wonder what is so special about these particular processes, that makes it possible to find 
the Wiener-Hopf factorization explicitly? It turns out that in all cases the Laplace exponent, defined as 
ip(z) = ln(E[exp(zXi)]), has some analytical structure which allows to factorize it as a product of two 
functions, which are analytic in the left/right complex halfplane. It is not surprising that the analytic 
structure of ip(z) plays such an important role, as there is a close connection between Wiener-Hopf 
factorization and the Riemann boundary value problems, see [12], [18] and the references therein. For 
example, if the process has hyper-exponential jumps [7], then ip(z) is a rational function and if X is a 
meromorphic process then ip(z) is a meromorphic function of a very special type, in both cases these 
functions can be easily factorized as products of two functions. One can formulate a general "meta- 
theorem": Wiener-Hopf factorization can be obtained explicitly if and only if ip(z) can be extended to a 
meromorphic function in the left or right complex halfplane. This principle helps to explain why no one 
has yet produced an explicit Wiener-Hopf factorization for one of the processes which are widely used in 
mathematical finance, such as VG, CGMY/KoBoL or generalized tempered stable processes (see [9] and 
the references therein for more information about these families of Levy processes). It turns out that in 
all these cases the Laplace exponent has a logarithmic or algebraic branch point in the complex plane, 
and, therefore, cannot be extended meromorphically. At the same time, we can use this meta-theorem 
to produce a large class of processes for which there is some hope to have an explicit Wiener-Hopf 
factorization: if the process has bounded jumps then it follows quite easily from the Levy-Khintchine 
formula that the Laplace exponent ip(z) is an entire function, and it might be possible to factorize it as 
a product of two functions and obtain some useful information about the Wiener-Hopf factorization. 

In this paper we consider a more general class: Levy processes with bounded positive jumps. There 
are two main reasons, one theoretical and one more practical, why we are interested in studying this 
class of processes. First of all, one can see that this is a very large class. In a certain sense it is "dense" 
in the class of all Levy processes: clearly, any Levy measure can be approximated arbitrarily close by 
truncating it at a large positive number. Therefore studying the Wiener-Hopf factorization for this class 
will lead to a better understanding of related results for general Levy processes. The second reason is that 
there are several situations where processes with bounded positive or negative jumps would be natural 
candidates for modeling purposes. One important example is ruin problem for the insurance company 
which is protected by the reinsurance agreement. In this case the size of each claim is essentially capped 
at a fixed level, and the amount of the claim above this level is being covered by the reinsurer. The 
value of the insurance company can be conveniently modeled by a spectrally negative Levy process with 
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bounded jumps, and now we have an interesting problem of how to compute numerically such important 
quantities as the ruin probability, discounted penalty function, etc. 

It is instructive to draw a parallel with the results of Lewis and Mordecki [29] on processes with 
positive jumps of rational transform (see also recent paper by Fourati [12] on double-sided exit problem 
for this class of processes). In their case the Laplace exponent of the ascending ladder process K,(q, z) 
(see [25] for the definition of this object) is a rational function, with all singularities in the left half-plane 
Re(z) < 0. In our case it will turn out that rz(q, z) is an entire function of a very special type: it belongs 
to the so-called Cartwright class (see [27] and the proof of Theorem 1). This makes it possible to factor it 
as an infinite product and to identify the Wiener-Hopf factors. There are also some similarities between 
the analytical structure for Levy processes with bounded positive jumps and meromorphic processes [21]. 
In both cases the positive Wiener-Hopf factor is given as an infinite product involving the solutions to the 
equation ip(z) = q in the half-plane Re(z) > 0. The major difference is that in the case of meromorphic 
processes the solutions to the equation ip(z) = q are all real and simple, while they are complex when 
the process has bounded positive jumps, and this fact makes the analytical theory more interesting and 
the computations somewhat more challenging. 

The paper is organized as follows: in Section 2 we present our main results on the Wiener-Hopf 
factorization for processes with bounded positive jumps, and we obtain an expression for the Wiener- 
Hopf factors as an infinite product in terms of the solutions to i/j(z) = q. We also study the asymptotics of 
these solutions, which will turn out to be very important for applications and numerical computations. 
In Section 3 we consider the spectrally negative case, and we obtain a series representation for the 
scale function W^ q '(x). A brief discussion of numerical methods and the results of several numerical 
experiments are presented in Section 4, while Section 5 contains the proofs of all results. 



2 Processes with bounded positive jumps 

Let us first introduce some notations and definitions. The Levy measure of the process X will be 
denoted by n(dx), and we will use the following notations for its tails: n + (a;) = n((a;, oo)) for x > 
and n~(x) = n((— oo,x)) for x < 0. In this paper we consider the class of processes with bounded 
positive jumps, thus we will assume that the Levy measure n has support on (— 00,/c]. Here k is the 
right boundary of the support of n, that is 

k = ini{x > : fl + (x) = 0}. (1) 

We will also assume that k > 0, so that we exclude the spectrally negative case, which will be considered 
in the next section. Note that at this stage we do not impose any restrictions on the Levy measure on 
the negative half-line. 

The Laplace exponent of the process X is defined as ip(iz) = ln(E[exp(izXi)]) for zGl, and it can 
be expressed by the Levy-Khintchine formula as follows 

k 

^(z) = l -a 2 z 2 + ^z+ J (e zx - 1 -zxh{x))U(dx), (2) 

—00 

where a > 0, /i G R and h(x) is the cutoff function. When the process has jumps of bounded variation, 
or equivalently, when 

1 

J |x|n(dx) < 00, (3) 
-1 
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we will assume that h(x) = 0, then \i can be interpreted as the linear drift of the process. When the 
jump part of the process has infinite variation, or equivalently, when condition (3) is violated, we will 
assume that h(x) = lr x> ^iy (or h(x) = 1 if E[|Xl|] is finite). Note that formula (2) implies that the 
Laplace exponent ip(z) can be analytically continued into the half-plane Re(z) > 0. Also note that ip(z) 
is real when z > 0, and that ip(z) = ip{z). In particular, the last property implies that if for some q G K 
the number z G C is a solution to the equation ip(z) = q, then so is z. 

Everywhere in this paper we will denote the first quadrant of the complex plane as 

Qi := {z G C : Re(z) > 0, Im(z) > 0}, 

and we will always use the principal branch of the logarithm and the power function, that is the branch 
cut will be taken along the negative half-line and for all z G C we have arg(z) G (— tt, 7r]. 



2.1 Analytic properties of the Wiener-Hopf factors 

The following theorem is our first main result. It describes the analytic structure of the Wiener-Hopf 
factors for processes with bounded positive jumps. 

Theorem 1. Assume that q > 0. Equation ip{z) = q has a unique positive solution Co and infinitely 
many solutions in Qi, which we denote by {Cn}n>i- Assume that ( n are arranged in the order of increase 
of the absolute value. The following statements are true: 

(i) Co has multiplicity one and Re(Cn) > Co for all n > 1. 

(ii) The series ^2 n>1 Re converges. 

(Hi) All of the numbers {( n }n>i, except possibly those of a set of zero density, lie inside arbitrary small 
angle it/2 — e < arg(z) < ti/2, and the density of zeros inside this angle is equal to 

lim #{C» : \Cn\ < r and tt/2 - e < arg(C n ) < tt/2} = k_ 



(iv) The Wiener-Hopf factors can be identified as follows: for Re (2) > 

' 0+(i.) := E [ e -**. W ] = e ¥ (l + ]J (l + ^ 

i-(-iz) := E [e zI -M] 



2tt 



z 

1 + S 



-1 



(5) 



q-ijj{z) <f>+(-iz)' 



The proof of Theorem 1 can be found in Section 5. 



Let us present a very simple example, which will illustrate the results presented in Theorem 1. 
Consider a process X t = kN t , where N t is the standard Poisson process. It is clear that X is a process 
with bounded positive jumps, and that its Laplace exponent is ip(z) = exp(kz) — 1. Solving equation 
ip{z) = q for q > we find that 

_ ln(l + q) _ ln(l + q) 2nvri 
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It is an easy exercise to verify that the series ^2 n>l Re (£~ converges, thus we have checked part (ii) 
of the Theorem 1. Next, all the zeros belong to the vertical line Re(z) = (o, they are equidistant and 
the spacing between them is equal to 2ir/k. This confirms statement (iii): all zeros (except for a finite 
number) lie inside arbitrary small angle 7r/2 — e < a,rg(z) < 7r/2, and the density of zeros inside this angle, 
which is inversely proportional to the spacing, is equal to k/(2ir). Finally, since X is a subordinator, we 
have S'e(g) = X e ( q ), thus 

r c -i -n-i r v'y i Q Q k * sinh (I ln(l + 5) ) 

E \e~ s °m] = E Te" c («)l = = = e^" — — — 

L J L J q-ip(-z) l + q-e- kz sinh (\{kz + ln(l + <?))) ' 

and we see that the infinite product representation in (5) is equivalent to the well-known infinite product 
formula for the hyperbolic sine function. 

It is also easy to verify the validity of Theorem 1 for a more general class of processes with double- 
sided jumps. Let us assume that for some h > the measure Il(dx) is supported on a finite subset 
of a lattice hZ, that is there exist m,l E N such that the support of II(dx) is equal to {—mh, — (m — 
l)h, . . . , —h, h, ...,(/ — l)h, lh}. In this case the right boundary of the support of the Levy measure 
is k = lh. Let X be a compound Poisson process defined by the measure II (da;) (note that X can be 
constructed as a linear combination of m + I independent Poisson processes). From the Levy-Khintchine 
formula (2) we find that the Laplace exponent of X is given by 

m I 

V>(z) = $>({-#}) [eT ihz - 1) + j>(07»}) {e? hz - 1) . 
3=1 i=i 

Note that the function ip(ln(w)/h) is a rational function, therefore using the change of variables z = 
ln(w)/h the equation ip{z) = q can be transformed into a polynomial equation of degree m + I. It is 
possible to prove (we leave it as an exercise) that this polynomial equation will have exactly m solutions 
inside the open unit circle {w G C : \w\ < 1} and exactly I solutions {u>i,. . . ,wi} outside the closed 
unit circle. The solutions ( n to the original equation if>(z) = q can now be found by solving equation 
exp(hz) = Wj, thus they are given by 



+ ; nez, i< 3 <i 



h h 

Again, it is easy to check that the series ^ n>1 Re (C^ 1 ) converges. Also, the solutions lie on I vertical lines, 
therefore all of them (except for a finite number) lie inside arbitrary small angle ir/2 — e < arg(z) < 7r/2, 
and the density of zeros inside this angle is equal to I x h/{2^) = k/(2ir). Infinite product representation 
for the positive Wiener-Hopf factor (5) is again equivalent to elementary infinite product expressions for 
certain trigonometric functions. 

In the above two examples we were able to describe the solutions to the equation ip(z) = q in a 
very precise form, but in the general case this will be a transcendental equation and there is little hope 
to obtain as much information about ( n . However, as our next result shows, we can obtain some very 
useful information about the asymptotic behavior of ( n , provided that the Laplace exponent ip{z) has 
regular growth as z — > 00. The connection between the regularity of growth and distribution of the 
zeros of an entire function is well-known, see, for example, chapters 2 and 3 in [27]. The main idea is 
that analytic functions which grow regularly at infinity also enjoy certain regularity in the distribution 
of zeros (and in fact the opposite is also true). In the following Theorem we impose a rather strong 
regularity condition on the growth of the Laplace exponent in the half-plane Re(z) > in order to 
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obtain an explicit asymptotic approximation for the solutions to ip(z) = q. This asymptotic expression 
for Cn would prove to be very useful in the next Section, when we will derive a series representation for 
the scale function W^ q '{x), and later in Section 4, when we will discuss numerical algorithms. 



Theorem 2. Assume that 



ip{z) = Ae kz z~ a + Bz b + o (e kz z- a ) + o (z b ) 



(6) 



as z — > oo in the domain Q\. Let us also assume that a > and b > 0. Then all sufficiently large 
solutions to ip(z) = q are simple and there exists m G Z such that 



n+m 



+ 



hi 



arg 



B 
A 
B 
A 



a + b) In 
1 



2nir 
~k~ 



(7) 



a + b) + 2n + 1 



7T 



0(1) 



as n — > +oo. 

The proof of Theorem 2 is presented in Section 5. 

Remark 1. Note that formula (7) implies that Q ri = (a + b) ln(n)//c + 2^n\jk + 0(1) as n — > oo, which 
again confirms statements (ii) and (iii) of Theorem 1: the zeros cluster "close" to the imaginary axis, or 
to say it more precisely, arg(£ ra ) n/2 as n — > oo; and secondly, the density of zeros in Q\ (which is 
inversely proportional to the average spacing between them) is equal to k/(2ir). 

Condition b > in Theorem 2 implies that ipiiz) -> oo as z oo, z e M, therefore X cannot be 
a compound Poisson process (see Proposition 2 in [3]). This shows that the two examples considered 
on page 4 do not satisfy the conditions of Theorem 2, but if we take these compound Poisson processes 
and add a drift (or Brownian motion with drift) then it is easy to check that the Laplace exponent of 
this perturbed process will satisfy (6). A natural question then is to describe sufficient conditions on 
the triple {fi,a,H}, which defines the Laplace exponent via (2), which will ensure that ip{z) satisfies 
asymptotic relation (6). Below we present a set of sufficient conditions. 

Definition 1. We will say that a real function f(x) is piecewise n-times continuously differentiable on 
an interval [a, b] if there exists a finite set of numbers {xk}i<k<m, such that 

(i) a = xi < x 2 < ■ ■ ■ < x m = b, 

(ii) f e C n ([a,b]\{xi,x 2 , ■ ■ ■ ,x m }), 

(iii) for each j — 0, 1, . . . , n and k = 1, 2, . . . , m 



there exist left and right limits f^\xk~) and f^\xk+) 



We will use the notation / G VC n [a, b}. In the case of an open interval (a, b) the definition of VC n (a, b) 
is very similar, except for condition a < x\ and x m < b. Similarly, one can define the remaining cases of 
intervals (a, b] and [a, b). 



Definition 2. We say that a Levy measure is regular if the following two conditions are satisfied: 
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(1) There exist constants C,a and {Cj, &j}i<j<rh such that 

m 

fi-(x)-C\x\- & -J2Cj\x\- & * = 0(1), z-KT, (8) 
where a, dj G (—00, 2) \ {0, 1} and acj < a. 

(2) There exists 116NU {0} such that 

(2a) for some constants C, a and {Cj,atj}i<j< m we have 

m 

n+ (z) - Cx- a - c 3%~ aj e PC n+1 [0, (9) 

where a, <x,- G (—00, 2) \ {0, 1} and aj < a; 
(2b) n+( n )(A;-) ^ 0; 

(2c) n + (x) G C n_1 (IR + ) (this condition is not needed for n = 0). 
Remark 2. Note that conditions (1) and (2a) imply that the Blumenthal-Getoor index 

0(11) = inf 1 7 > : J \x\m(dx) < 00 j (10) 

is equal to /3(II) = max(a, a, 0). 

Definition 2 is not very easy to interpret, therefore we will try to give some intuition behind these 
conditions. Conditions (1) and (2a) guarantee that the Levy measure is sufficiently well-behaved in the 
neighborhood of zero. This will help us to ensure that the main term of if)(z) grows as z — > 00 exactly as 
z b , and does not contain any logarithmic terms. Conditions (2b) and (2c) are slightly harder to interpret. 
Essentially, they imply that the Levy measure restricted to M + has its "worst" possible singularity at 
the right-end point of its support. Let us consider the following example, where conditions (2b) and (2c) 
are violated. 

Example 1. Assume that the Levy measure is given by 

U(dx) = l {x<0 }e x dx + l{ 0<x< 4}dx + S 3 (dx), (11) 

so that n(dx) has an atom of mass one at x — 3. Because of the atom at x = 3 we know that Yl + (x) is 
not continuous, therefore we are forced to take n = in the Definition 2. But since we have no atom at 
x = k = 4, we find that Tl + (k— ) = 0, which violates condition (2b), thus we conclude that the measure 
n(dx) is not regular. 

Next, let X be a process which has a Levy measure (11) and linear drift \x = 1. We will check that 
the Laplace exponent of the process X does not satisfy (6). We compute the Laplace exponent using 
the Levy-Khintchine formula (2) and find that it has asymptotics 

4>(z) = z+ — + e 3z + 0(l), 
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as z — > oo, Ke(z) > 0. Now, it is easy to see that in the domain < Ke(z) < ^ In | ^ | we will have 
ip(z) = z + e 3z + o(e 3z ), while in the domain Re(z) > 21n|z| we'll have ip(z) = e 4z /z + o(e Az /z). 
This implies that we cannot find a single uniform asymptotic formula for tp(z) as in (6). This happens 
because the "worst" singularity of IT(da;), which is the atom at x — 3, is not located at the right boundary 
x = k = 4. One can also check that asymptotic expression (6) will be satisfied if we replace S^(dx) by 
S^(dx) in (11) or if we add a second atom at x = 4, and at the same time this will also give us a regular 
Levy measure according to the Definition 2. 

In the following example we exhibit a large family of regular Levy measures (it is an easy exercise to 
verify all the conditions of Definition 2). 

Example 2. Assume that the Levy measure Il(dx) has a density tc(x) given by 

tt(x) = l^oj/CaOlxr 1 - 6 + l {0<a ,<fc}/(a:)a:~ 1 ~ a , (12) 

where a, a G (—oo, 2) \ {0, 1} and functions /, / satisfy the following conditions: (i) f(x) and f(x) can 
be represented by convergent Taylor series in some neighborhood of zero; (ii) f(x) is VC [0,A;]; (iii) 
f(k-) > 0. Then n(dx) is regular. 

The above example shows that there are indeed many interesting Levy processes with regular Levy 
measure. For example, we can take one of the widely used processes in mathematical finance, such 
as CGMY/KoBoL or generalized tempered stable (see [9]), truncate its Levy measure at any positive 
number and we will obtain a regular Levy measure. 

The next Proposition shows that if the Levy process has a regular Levy measure, then its Laplace 
exponent satisfies the asymptotic expansion (6) and therefore the roots ( n have simple asymptotic ap- 
proximation given by (7). 

Proposition 1. Assume that X is not a compound Poisson process and that the Levy measure of X is 
regular. Then the asymptotic expression (6) is true, with parameters A = (— l) n fl + ( n )(&;— ), and a = n. 
The parameters B and b can be identified as follows: 

(i) if a > 0, then B = cr 2 /2 and b = 2, 

(ii) if the process has paths of bounded variation and /x ^ 0, then B = fi and 6=1. 

In the remaining cases, when the process has paths of unbounded variation and a = 0, or when the process 
has paths of bounded variation and fi = 0, we have 6 = /3(H) = max(a, a) and 

{-Ce- wia T(l- a), if a > a, 

-CT(l-d), if a < a, 

-(Ce~™ + C)T(1 - a), if a = a. 

Moreover, the asymptotic expression for ip' (z) can be obtained from (6) by taking derivative of the right- 
hand side. 

The proof of Proposition 1 can be found in Section 5. 
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2.2 Partial fraction decomposition and distribution of S e ^ 

As we have mentioned in Section 1, in the case when X has positive jumps of rational transform the 
positive Wiener- Hopf factor is a rational function (see [29]). By performing the partial fraction decom- 
position of this rational function and inverting the Laplace transform one can obtain the distribution of 
S e ( q ) explicitly. The same procedure works for meromorphic processes (see [21]), though in this case we 
must work with meromorphic functions instead of rational functions, and things become slightly more 
technical. In our case, when the process has bounded positive jumps, formula (5) tells us that the pos- 
itive Wiener- Hopf factor <pg(iz) is also a meromorphic function, thus we might hope to follow the same 
procedure and obtain a series representation for the distribution of S e ( q ), which would be very useful for 
applications. Unfortunately it turns out that proving existence of the partial fraction decomposition for 
4>q(iz) of the form (5) is a much harder problem, and we were not able to give a completely rigorous 
proof of such a result or to find such a result in the existing literature. However, if one assumes that 
such a partial fraction decomposition exists, it is rather easy to obtain the form of its coefficients and the 
resulting expression for the distribution of S e f q y Let us sketch here the main steps, and later, in Section 
4 we perform several numerical experiments which seem to confirm our conjecture. 

First of all, it is very likely that all the zeros {Cn}n>i of the function if)(z) — q are simple. This 
was proved in Theorem 2 for large zeros and for Co- For other roots one could use the following (rather 
informal) argument. We know that z is a solution of ip(z) = q with multiplicity greater than one if and 
only if i/j'(z) = 0. We can rephrase this statement: equation ip(z) = q has a solution of multiplicity 
greater than one if only if q = where C is a root of ip'(z). But ip'(z) is analytic in the half-plane 

Ke(z) > 0, thus it has a discrete set of zeros, and for any given complex root of ip'(z) it is very unlikely 
that VKC) wm De a re& l positive number, thus it is very unlikely that there exists q > such that tp(z) — q 
has a multiple root. Even if such a value of q exists, we see that in the worst possible case, there can be 
only finitely many such values on any compact subset of R, this shows that the assumption that all the 
zeros {Cn}«>i ar e simple should not be very restrictive for practical purposes. However, of course this 
fact would require a rigorous proof in future work. 

Next, assuming that is regular for (0, oo), or equivalently, that the distribution of S e ^ has no atom 
at zero, we must have 4>q(iz) = E[exp(— zS e ^)] — > as Re(z) — > +oo, thus it is reasonable to expect 
that the partial fraction decomposition for <pti}z) should be of the form 



(iz) = E [e~ zS -M] 



k z 

e 2 



1 + 



Cc 



n>l 



n i+ 



1 + 



Cn 



a 



n>l 



+ 



_Z + Cn Z + Cn. 



(13) 



where a n = Res(0+ (iz) : z = — Cn) f° r ^ > 0. Using the above infinite product representation we can 
easily compute the residues at points ( n and obtain 



a = Coe~^ Co J] 



m>l 



Co 



and 



Co|U 2 e 



2 p o 



2Im(C„)(C„ - Co 



n 

m>l 



Cn 

Cm 



Cn 



Therefore, provided that there exists a partial fraction decomposition of the form (13), we can use the 
uniqueness of Laplace transform and conclude that the density of the supremum S e ( g ) should be given 
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by the following infinite series 

^(S c{q) <x) = a e-^ x + 2 ^ Re [a n e~^ x ] . (14) 

n>l 

Again, we emphasize that at this point the existence of a partial fraction decomposition (13) is just a 
conjecture, which would have to be proven rigorously in future work. 



3 Scale functions for spectrally negative processes with bounded 
jumps 

Everywhere in this section we will assume that Y is a spectrally negative Levy process with bounded 
jumps, so that the Levy measure Hy(da;) is supported on the interval [—k, 0) where k > 0, and k is the 
smallest such number. From the Levy-Khintchine formula (2) it follows that in this case the Laplace 
exponent ipy(z) = lnE[exp(zVi)] is given by 



u 

iPy{z) = + l^z + J (e zx - 1 - zxh(x)) Uy(dx) 



and we see that ipy(z) is an entire function which is convex for real values of z. Since ipy(0) = 0, it is 
clear that for each q > the equation i/jy(z) = q has a unique positive solution z = $((?), and in fact 
it is known from the general theory of spectrally negative processes that this is a unique solution in the 
half-plane Re(z) > 0, see chapter 8 in [25]. 

For q > the scale function W ( - q \x) is defined as follows: W' 9 '(x) = for x < and on [0, oo) it is 
characterized via the Laplace transform identity 



/ 



e- zx W ig \x)dx = — — , Re(z) > $(q). (15) 

ip Y {z) -q 



The scale function can be considered as the main building block for the vast majority of fluctuation 
identities for spectrally negative processes, see [25, 23] for many examples of such identities. Here we 
will present one fundamental identity, which is related to the exit of the process Y from an interval. If 
we define the first passage time = inf{t > : Y t > a} and similarly Tq~ = inf{t > : Y t < 0}, then 
Theorem 8.1 in [25] tells us that 



E 



X 



6 i K + <r -} 



W^Xx) 



In fact, this identity justifies the name "scale function": we see that W^ q \x) plays an analogous role to 
scale function for diffusions. 

Our main goal in this section is to obtain an expression for the scale function W^ q \x) in terms of the 
Laplace exponent and the roots of the entire function ipy{z) — q. We will consider spectrally negative 
Levy processes, whose Levy measure satisfies the following definition. 

Definition 3. We say that the Levy measure of a spectrally negative process Y is regular if there exists 
nGNU {0} such that 
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(a) for some constants C, a and {Cj, cij}\<j<m we have 

m 

fl-(x) -C\x\- a -J2 C i\ x \~ ai e VC n+1 [~k,0}, 

3=1 

where a, atj G (—00, 1) U (1, 2) and a, < a; 

(b) n-(")(-A; + ) ^ 0; 

(c) n~(x) G C n_1 (M _ ) (this condition is not needed for n = 0). 

By considering the dual process K = — K and using Proposition 1 we see that if Y has a regular Levy 
measure then its Laplace exponent satisfies 

ih{-z) = Ae kz z- a + Bz b + o (e kz z- a ) + o (z b ) (16) 

as z — > 00 in the domain Qi, where a > and b > 0. Moreover, the parameters in the asymptotic 
expression (16) are given by 



( A = U- ( - n H-k + ), and a 



n, 



£ = and 6 = 2, if a > 0, 

2 (17) 
5 = — /i, and 6 = 1, if u = and a < 1, 

£ = -Ce -,ria r(l - a), and 6 = a, if <r = and a > I. 

Theorems 1 and 2 tell us that equation ipy(— z) = q has infinitely many solutions {( n }n>o in Qi, 
moreover the first solution £0 is rea l an d positive, we have ( n > Co fo r n > 1 and the large solutions £ n 
satisfy asymptotic relation (7) with constants A, a, B and b as given in (17). 

The next Theorem is our main result in this section. It provides a series representation for the scale 
function W ( - g \x) in terms of the Laplace exponent ipy(z) and the numbers ( n . Its proof can be found in 
Section 5. 

Theorem 3. Assume that q > or q = and <$>(q) > 0. If the Levy measure of Y is regular and all 
solutions to ipy{z) = q are simple, then for x > we have 



^Rn)J ' 



Wi9) ( x ) = T, ^ \\ + 777 T\ + 2 E Re 
^($(g)) Vr(-Co) ^ 

where the series converges uniformly on [e, 00) for every e > 0. 

Formula (18) is in fact very similar to the corresponding expression for the scale function for mero- 
morphic processes, see [24] and [23]. 

In the very unlikely case that some of the solutions to ^>y{z) = q are not simple (see the discussion 
on page 9) the expression in the right-hand side of (18) would have to be modified. The coefficients 
exp(— Cnx)/ip' Y (~ ( n ) are the residues of exp (zx) /(tpy(z) — q) at z = —( n (see the proof of Theorem 3 in 
Section 5), and these coefficients would have to be appropriately modified if ( n is a root of ipy(— z) — q 
of multiplicity greater than one. 
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4 Numerical examples 



The main reason why we are interested in the analytical structure of the Wiener-Hopf factorization is 
that its understanding can lead to efficient numerical algorithms for computing such important objects 
as the distribution of supremum S e ( q ) or infimum I e (q), or the scale function W^ q \x). These are not easy 
problems for general Levy processes. For example, computing the distribution of S e r q ) in general involves 
evaluating numerically two integral transforms: first one has to compute the positive Wiener-Hopf factor 
via the formula (see [29]) 



E \e izS ^] = exp [— / In ( *—) d " 

1 J |_27ri Jm \q-^(iu)J 



q — ip{iu) J u(u — z) 



lm(z) > 0, 



and then perform an inverse Fourier transform to recover the distribution of S e ( q y Similarly, computing 
the scale function W^ q '(x) in general is equivalent to inverting Laplace transform in (15) (see [23] for the 
detailed discussion and comparison of several numerical algorithms for computing the scale function). 

The results presented in this paper lead to quite different approach for computing the distribution of 
S e ( q ) and the scale function W^ q '(x). This approach does not rely on the numerical evaluation of multiple 
integral transforms; instead, the main ingredients are the solutions to the equation ip(z) = q and infinite 
series representations (14) and (18). Therefore, this approach is very close in spirit to the techniques 
that are used for processes with positive jumps of rational transform [29] or meromorphic [21] processes. 

Our main goal in this section is to give a brief description of the numerical algorithms and techniques 
which are suitable for Levy processes with bounded positive jumps, in particular, we would like to show 
that series expansions (14) and (18) may lead to efficient numerical computations. However, the detailed 
investigation of these numerical algorithms is beyond the scope of the current paper and we will leave to 
future work such important questions as the speed of convergence, rigorous error analysis, etc. 



4.1 Preliminaries 

In order to implement the numerical algorithms based on formulas (14) and (18) we have to solve the 
following problem: how can we find the complex solutions of the equation ip{z) = g? As we will see, it 
is not an easy problem, yet it can be solved rather efficiently provided that we use the right techniques. 

The main problem in finding the zeros of ip{z) — q is that all of them (except for ( ) are complex 
numbers. Note that the real zero (o can be easily found by bisection method followed by Newton's 
method. The large zeros of ip(z) = q satisfy asymptotic relation (7), thus they can be found by Newton's 
method which is started from the value in the right-hand side of (7). The only problem that remains is 
how to compute the complex zeros of ip(z) — q which are not too large. 

The problem of computing the zeros of an analytic function inside a bounded domain has been 
investigated by many authors, see [10, 34, 35] and the references therein. We will follow the method 
presented in [10], which is based on Cauchy's Argument Principle. Let us recall this important result. 
We denote the change in the argument of an analytic function f(z) over a piecewise smooth curve C as 



Aarg(/, C) = Im 



c 



m z . 



(19) 



provided that f(z) ^ for z G C. Cauchy's Argument Principle states that if C is a simple closed 
contour which is oriented counter-clockwise, and f(z) is analytic and non-zero on C and analytic inside 
C, then Aarg(/, C) = 2ttN, where N is the number of zeros of f(z) inside contour C. 
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This result leads to a practical iterative procedure to determine the complex zeros of an analytic 
function f(z) inside a closed contour. We start with an initial rectangle R in the complex plane and 
proceed through the following sequence of steps: (i) compute N = N(R) - the number of zeros of f(z) 
inside rectangle R, (ii) if N — 0, then we stop, (iii) if N — 1 - we try to find the zero using Newton's 
method started from a point inside the rectangle, (iv) if N > 1 or if the Newton's method in step (iii) 
fails - then we subdivide the rectangle R into a finite number of disjoint rectangles Rj. For each of 
the smaller rectangles Rj we proceed through the same sequence of steps (i)— >(iii)— >(iv). At some 
point the rectangles which contain zeros of f(z) become very small, therefore the starting point of the 
Newton's method is close to the target and the Newton's method converges to the zero of f(z). This 
shows that in theory we should be able to recover all zeros of f(z) inside R. There are some technical 
issues which arise when at some step of the algorithm we obtain a rectangle with one of the zeros being 
extremely close to the boundary of this rectangle, however we found that in our numerical experiments 
this issue was not a big problem. More details of this algorithm and some numerical examples can be 
found in [10]. 

Next, let us review some basic facts about the incomplete gamma function, which will be used 
extensively in this section. The incomplete gamma function is defined for Re(s) > and z > as 

7M = /«-V«*, (20) 



It is known (see section 8.35 in [15]) that the function z i— >■ z~ s, -f(s, z) is an entire function and it can be 
represented by a Taylor series (which converges everywhere in the complex plane) as follows 

Here, as usual, (a) n = a(a + l) . . . (a+n — 1) denotes the Pochhammer symbol and the function i-F\(a, b; z) 
which appears in the above equation is the confluent hypergeometric function 

^M-£||S, (22) 



n>0 



see chapter 6 in [11] for an extensive collection of results and formulas related to this function. 

While the incomplete gamma function is not one of the elementary functions, it can still be easily 
evaluated everywhere in the complex plane. In fact, numerical routines for evaluating this function are 
provided in such computational software programs as Maple and Mathematica. One should use different 
strategies for computing 7(5,2) depending on whether \s\ and/or \z\ is large. However, we will only need 
to compute j(s, z) for a fixed s (which is real and not large) and for various values of z. In this case 
the numerical algorithms are based on one of the infinite series expansions presented (21) when \z\ is 
not large, or on the various asymptotic approximations (and expansions in continued fractions, such as 
formula 8.358 in [15]) when \z\ is large. We refer to [16] or [33] for all the details. 

Finally, we would like to mention that the code for all numerical experiments was written in C++ 
and the computations were performed on a standard laptop (Intel Core i5 2.6 GHz processor and 4 GB 
of RAM). 
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(a) (b) 

Figure 1: On graph (a) we present the first fifty roots ( n (crosses) and their approximations (circles) 
given by Theorem 2. On graph (b) we present the distance from ( n to its approximation (note that we 
use logarithmic scale for the y-axis). 



4.2 Numerical example 1: processes with double-sided jumps 

For our first numerical experiment, we consider a generalized tempered stable process (see [9]), also 
known as KoBoL process, with the Levy measure truncated at a positive number k 

Il(dx) = l {x<Q} Cae^ x \x\- l - & dx + l {Q<x<k] Cae- px x- l - a dx. (23) 

Proposition 2. Let X be a Levy process, with the Levy measure il(dx) given by (23). Then the Laplace 
exponent of X is given by 

if>(z) = -a 2 z 2 + fiz- CT(1 -&)0 + z f + Ca{(3 - z) a j(-a, k((3 - z)) + rj, (24) 

where r\ is chosen so that ^(0) = 0. 

The proof of Proposition 2 can be found in Section 5. Note that when z < (5 and k — > oo, then 
7(— a, k((3 — z)) — > r(— a) (this follows from (20) when a < 0, see also formulas (8.356.3) and (8.357.1) 
in [15]). This confirms the intuitively obvious result that as the cutoff k becomes very large, the Laplace 
exponent of the truncated process converges to the Laplace exponent of the generalized tempered stable 
process, see Proposition 4.2 in [9]. Note also that the function (/3 — z) a ^{— a, k(/3 — z)) which appears in 
(24) is an entire function of z, which confirms the fact that ip(z) is analytic in the half-plane Re(z) > 0. 
As in the case of generalized tempered stable processes, when cr > or q > 1 or d > 1 we have a process 
of infinite variation. Finally, in the case when o = and a < 1 and a < 1 we have a process of finite 
variation, and in this case fi corresponds to the linear drift of this process. 

We consider the following parameters: 

a = \ 1 /i = -2, (7 = 6=1, a = d = 0.5, = 1, f3 = 2, k = 1. 

These parameters give us a process with negative linear drift and finite variation infinite activity jumps. 
Note that since cr = 1 we have a process with paths of unbounded variation. We also set q = 1. 

First we compute 1000 roots ( n using the method discussed in Section 4.1. Overall, it takes just 
0.15 seconds to compute 1000 roots. The results are presented on Figure 1. Note that, as expected, the 
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Figure 2: The density of S e ^ for q — 1. 



approximation to £ n provided by (7) becomes better and better as n increases, but is not so good for 
small values of n. 

Next, we compute the density p(x) of the supremum S e ^ using the series representation (14). The 
results are presented on Figure 2. After we have pre-computed and stored the roots ( n , computing 2000 
values of p(x) takes just 0.26 seconds. In order to test the accuracy we have numerically computed the 
integral of p(x) on the interval [0, 10] (which should be close to P(S e (g) > 0) = 1). The result is equal 
0.985 if we use 1000 roots and 0.995 if we use 5000 roots. 

4.3 Numerical example 2: a family of spectrally negative processes 

For our second numerical experiment we consider a spectrally negative Levy process Y, with the Levy 
measure defined as follows 

LV(dx) = l^ k<x<0 }Cae ,3x \x\~ 1 ~ a dx. 

Note that the dual process Y = —Y belongs to the class of processes described in Proposition 2, therefore 
the Laplace exponent of Y is given by 

4> Y (z) = ^a 2 z 2 + fiz + Ca{(3 + z) a 7(~a, Hf3 + z)) + 77, (25) 

where again rj is chosen so that ^y(0) = 0. 
We fix the following values of parameters 

ae{0,l}, /i = 2, C = l, a = 0.5, (3 = 1, k = 1, (26) 

which define a spectrally-negative process with infinite activity /finite variation jumps and paths of fi- 
nite/infinite variation depending on whether a = or o = 1. We compute 1000 numbers {( n } using the 
algorithm presented in Section 4.1. This computation takes 0.06 seconds. The qualitative behavior of 
the roots ( n is very similar to the one presented on Figure 1. Computing the scale function via series 
representation (18) is also very fast: it takes just 0.07 seconds to compute 1000 values of W^ q \x) for 
equally spaced points x G [0, 3]. 
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(a) a = (b) a = 1 

Figure 3: The scale function VK^(x) for the parameter set (26) 

The results of computations are presented on Figure 3. From Lemma 8.6 in [25] we know that 
W^(0) = if the process Y has infinite variation, and W^ q ' = if the process has bounded variation 
(where /i is the linear drift). One can see from Figure 3 that our numerical results are in perfect 
agreement with the theoretical prediction. It would also be very interesting to compare the accuracy and 
performance of this algorithm for evaluating W^ q \x) with the methods used in [23], however we have 
decided to leave this for future work. 

5 Proofs 

Lemma 1. Let v(dx) be a finite positive measure such that v([a,+oo)) = for some a el. Define 
k = inf{a G K : v([a, +oo)) = 0}. Then for every e > there exists £ = £(e) > such that for all z > £ 
we have 



e (k-e)z 



< J e zx v(dx) < e {k+e)z . (27) 



Proof. Assume that e > 0. Since v((k, oo)) = we have for all z > 

e- (fc+e)2 y e zx u(dx) = J e z{x ~ k - e) v{dx) < e~ ez u(R), 

H {-°o,k] 

and the right-hand side in the above inequality goes to zero as z — > +oo. This proves the upper bound 
in (27). Similarly, 

e - {k - e)z J e zx u(dx) > e - (k - e)z J e zx u(dx) > ei z u([k - e/2, k}). 

R [k-e/2,k] 

According to our definition of k, the quantity u([k — e/2, k]) is strictly positive, thus the right-hand side 
in the above inequality goes to +oo as z — > +oo, which proves the lower bound in (27). □ 
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Proof of Theorem 1: Let us prove (i). Levy-Khintchine formula (2) implies that the function ip(z) 
is analytic in the half-plane Re(z) > and convex for z > 0. Since ^(0) = and ip(z) increases 
exponentially as z — > +oo we conclude that there exists a unique simple real solution to ip(z) = q, which 
we will denote Co- 

Let us prove that there are no other solutions in the vertical strip < Re(z) < Co- From the definition 
of the Laplace exponent we find 

gtReW*)) = |e [ e zXt ] | < E [ e Re{2)Xt ] = e ^(R e ( z )). 

This shows that for < Re(z) < Co we have Re(ip{z) — q) < ip(Re(z)) — q < 0, therefore ip(z) ^ q in this 
vertical strip and we have proved the first part of Theorem 1. 

Next, let us prove (ii), (iii) and (iv). Let us consider the ascending ladder process (L _1 ,if) and its 
Laplace exponent n(q, z). Similarly let k(q, z) denote the Laplace exponent of the descending ladder 
process see section 6.2 in [25] for the definition and properties of these objects. Let A(d£, dx) 

denote the Levy measure of the bivariate subordinator (L _1 ,iJ) (see section 6.3 in [25]). One can check 
that n(q, z) can be expressed in the following form 

oo 

K(q, z) = K(q, 0) + az - J (e~ zx - l) A ((?) (dx) (28) 

o 

where a > and 

POO 

A^(dx) = / e~ qt A(dt,dx). 
Jo 

Formula (28) implies that the function z (->■ n(q, z) is the Laplace exponent of a subordinator with the 
Levy measure Aw(drc) and drift a, which is killed at rate «(g, 0). Note that A^°\dx) is the Levy measure 
of the ascending ladder height process H. The jumps in the process H happen when the process X 
jumps over the past supremum, thus it is clear that if the jumps of X are bounded from above by k, 
then the same is true for the process H. Therefore A(°)((fc,oo)) = 0, and since for each Borel set B the 
quantity AS q \B) is decreasing in q we conclude that hS q \{k, oo)) = for all q > 0. Therefore we have 
proved that the support of the measure A^) (dx) lies inside the interval (0, k] . 

Let us define k = ini{x > : A {q \(x, oo)) = 0}. Note that k < k, since we have established already 
that the support of the measure A^ q \dx) lies inside the interval (0, k]. Using formula (28) and the fact 
that A^ q \dx) has finite support we conclude that n(q, z) is an entire function of z. 

Using the Wiener-Hopf factorization (see Theorem 6.16 in [25]) we find that for Re(z) < 

E [e zS ^i] | < E [ e Re ( z ) 5 cM] < 1, (29) 

this shows that the function n(q, —z) has no zeros in the half-plane Ke(z) < 0. By the same argument 
we conclude that k(q, z) has no zeros in the half-plane Ke(z) > 0. Therefore, using the Wiener-Hopf 
factorization q — ip(z) = rz(q, —z)k(q,z) we find that all zeros of n(q, —z) in the half-plane Re(z) > 
coincide with the zeros of q — ip(z). Recall that we have labeled these zeros as {Co, Cm Cn}n>i, where 
Cn £ Qi are arranged in the order of increase of absolute value. 

Next, we use Proposition 2 in [3] and the fact that n(q, z) is the Laplace exponent of a subordinator 
to conclude that n(q, iz) = O(z) as z — > oo, z 6 R. Therefore the following integral converges (which is 



ft(g,0) 
K(q, -z) 
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equivalent to saying that n(q, iz) belongs to Cartwright class of entire functions) 



max(ln(|/t(g,iz)|),0) 

1 + Z l 

R 

and we can apply Theorem 11, page 251 in [27] (see also remark 2, page 130 in [28]) to conclude that 
n(q, —z) can be factorized as 

: \ j-j- / _ : \ / , : 

n>l 



-fe-.)=«(*0),..(l- 5 J! l-j; l-j). (30) 



Moreover, all of the roots ( n (except for a set of zero density) lie inside arbitrarily small angle 7r/2 — e < 
arg(z) < 7r/2, the density of the roots in this angle exists and is equal to kj (2ir) and the series R- e (CT l ) 
converges. 

Using (28) and the following result 

E[e~ zS ^] =J ) taz)- K{q > 0) 



(see Theorem 6.16 in [25]) we obtain formula (5) for the Wiener-Hopf factors, and in order to finish 
the proof we only have to show that k = k. This fact seems to be intuitively clear, as it means that 
the upper boundary of the support of the Levy measure of the process X is exactly equal to the upper 
boundary of the support of the Levy measure of the ascending ladder height process H . However, we 
were not able to find a simple probabilistic argument to prove this statement, and we will use an analytic 
approach instead. Assume that k < k and define e = (k — k)/3. Using the Levy- Khint chine formula (2) 
and Lemma 1 one can check that IVK- 2 )! > exp((fc — e)z) for all z > large enough. Similarly, using (28), 
Lemma 1 and the fact that A'^(cLe) has support on (0, k] we can check that \n(q, —z)\ < exp((k + e)z) for 
all z > large enough. From the Wiener-Hopf factorization k(q, z) — (q — ip(z))/K(q, —z) we conclude 
that 

e (k-e)z 

for all z > large enough. But this is not possible, as we know that k(q, z) is the Laplace exponent 
of a subordinator, thus k(q,z) = O(z) as z — > oo, Re(z) > (see Proposition 2 in [3]). Therefore the 
inequality k < k is not true. At the same time we must have k < k, since the support of the measure 
A")(dx) lies inside the interval (0, k). This implies that k = k and ends the proof of parts (ii), (iii) and 
(iv) of Theorem 1 . □ 



Recall that Aarg(/, C), which was defined in (19), denotes the change in the argument of f(z) over 
a curve C. The following result will be used in the proof of Theorem 2. 

Proposition 3. Assume that f and g are analytic on a piecewise curve C and f(z) ^ for z G C . If 

for some e G (0, 1) we have \g(z)\ < e\f(z)\ for all z G C , then 

|Aarg(/ + s,C)-Aarg(/,C)| < 4e. (31) 
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Proof. From the definition of Aarg(/, C) (19) it follows that 

Aarg(/ + g,C) = Aarg(/, C) + Aarg(l + g/f, C). 

Due to the condition \g(z)/ f(z)\ < e for z G C we know that the set {u; = 1 + g(z)/ f(z) : z G C} lies 
inside the circle of radius e with center at one. Using elementary geometric considerations we check that 

max{|arg(w)| : \w — 1| < e} = arcsin(e) < 2e, 

thus the change of the argument of any curve lying inside the circle | w — 1 1 < e cannot be greater than 
4e, which proves (31). □ 



We will also need the following two facts: 

(i) For any real u > and any complex number v for which |arg(t>)| < 7r/2 we have 

|arg(w + t>)| < |arg(u)|. 

(ii) For any two complex numbers u, v we have 

\u + v | > cos(arg(w) — arg(t>))(|u| + \v\). 

Both of the above facts can be easily verified by elementary geometric considerations 



(32) 



(33) 



Proof of Theorem 2: Let us denote the set of solutions to if)(z) = q in Q 1 as Z = {Cn}n>i and 
introduce numbers 



k 



In 



(a + b) In 



2titc 



+ 



k 



m-g \ — J + (^(a + b) + 2n + 1 | rr 



First let us prove that every solution to the equation ip(z) = q which has sufficiently large absolute value 
must be close to one of z n . From the asymptotic expansion (6) we find that when z — > oo, z G Z we 
have 



.Mz 



B 
A 



z a+b (l + o(l)). 



Considering the absolute value of both sides of (34) we conclude 



Re(z) = - In 
k 



o ~\~ b .. 

H : — in \z\ + o(l), z — > oo, z £ Z. 

k 



From the above asymptotic expression it follows that Re(z) = 0(ln \z\), which in turn implies 

7T 

arg(z) = - + o(l), z ->• oo, z E Z. 
Next, considering the argument of both sides of (34), we find that 

arg (e fc2 ) = arg ( -— J + arg (z a+b ) + arg(l + o(l)) (mod 2vr), 



(34) 



(35) 



(36) 
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(a) (b) 
Figure 4: Illustrations to the proofs of Theorem 2 and Theorem 3. 



and using (36) we conclude that there exists an integer number n such that 



+ (2n + l)?r + -(a + b)n + o(l). 



klm(z) = arg 

Finally, asymptotic expressions (35) and (37) imply that 



(37) 



Re(z) 



1 

a- 



In 



+ (a + b) In 



2rt7r 



+ o(l) 



which together with (37) shows that every sufficiently large solution to ip(z) = q must be close to one of 
the numbers z n . 

Now let us prove the converse statement: for all n large enough there is always a solution of the 
equation ip(z) = q near a point z n . We set e = l/(16k(a + b)), assume that n is a large positive integer 
and consider the following contour L = L(n) = L\ U L 2 U L 3 U L 4 , defined as 



Li 


= In 






e C : 


Re(z 


u 


= L 2 




= 


G C : 


lm(z 


u 


= Ls 




= {* 


G C : 


Re(z 


u 


= L 4 




= {* 


G C : 


lm(z 



71 - 



7T 
7T 



71 - 
k 



As we see on Figure 4a, L is a rectangle of dimensions 27r/k and en, which contains exactly one point 
z n for n large enough. We assume that this contour is oriented counter-clockwise. Our goal is to prove 
that Aarg(?/>(z) — q,L(n)) = 2tt for all n large enough, and our strategy is to show that the change in 
the argument over L\, L 2 and L4 is small, while the change in the argument over L3 is close to 2tt. 
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First of all, it is clear that the number of zeros of ip(z) — q inside the contour L is the same as the 
number of zeros of F(z) = z a {ip(z) — q) inside the same contour. Asymptotic expression (6) tells us that 

F(z) = Ae kz + Bz a+b + o (e kz ) + o (z a+b ) , z oo, z G Qi. (38) 

Let us first consider the interval L\. Since arg(z) = ir/2 on this interval, we have Aarg(z a+b , L\) — 0. 
Equation (38) implies that F(z) = Bz a+b + o (z a+b ) when z G Li and n — > +oo, thus we use Proposition 
3 and conclude that for all n large enough we have |Aarg(F, Li(n))\ < 1/4. 

Let us consider the contour L 2 . From the definition of this contour it follows that for all z G L2 

arg^exp^-^(a + 6)jj =0. (39) 

Also, looking at Figure 4a one can check that n/2 — 6 < arg(z) < 7r/2 for all z G L2, where we have 
defined 



9 = 6(e, n) = arctan 



lm(z n ) -ir/kj ' 



Note that as n — > +00 we have 6(e,n) — > arctan(e/c/(27r)), and the latter quantity is smaller than ke. 
This implies that for all n large enough we have #(e, n) < ke. Thus we have proved that for all n large 
enough we have 7r/2 — ke < arg(z) < 7r/2 when z G L 2 , which is equivalent to 



arg ( z a+b ex.p ( -^(a + 6) 



2 

From (39), (40) and property (32) it follows that for all z G L2 the number 

w = — exp f kz - —(a + b) j + z a+b exp f -y(o + b) 
lies in the sector |arg(w)| < 1/16. From here we find that 



<k(a + b)e= zeL 2 . (40) 

lb 



\A&rg{Ae kz + Bz a+b ,L 2 )\ < I, (41) 



(42) 



and at the same time, with the help of property (33) we deduce that for all z G L 2 

\Ae kz + Bz a+b \ > cos (\Ae kz \ + \Bz a+b \) . 

Next, if g(z) = o (e kz ) + o (z a+b ) as z — >■ 00, then 

g(z) = o(\Ae kz \ + \Bz a+b \) , 

and we again can use (41) and (31) with f(z) := F(z) and g(z) defined above, to conclude that for all 
n large enough 

\Aarg(F(z),L 2 (n))-A&Tg(Ae kz + Bz a+b ,L 2 )\ < 1/8. 
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The above inequality and estimate (41) imply that for all n large enough we have |Aarg(F(2;), L,2(n))\ < 
1/4. Using exactly the same technique we obtain an identical estimate the change of argument over 
L 4 (n). 

Finally, on the contour L 3 we have F(z) = Aexp(kz) + o(exp(kz)). Since Aarg(exp(/cz), L 3 ) = 2tt, we 
use Proposition 3 and conclude that for all n large enough |A(F, L^{n)) — 2n\ < 1/4. Combining these 
four estimates we see that for all n large enough we have |A(F, L(n)) — 2tt\ < 1, and since we know that 
Aarg(F, L) must be an integer multiple of 2ir we conclude that Aarg(F, L) = 2tt, thus there is exactly 
one solution to ip{z) = q inside the contour L(n). From the first part of the proof we know that every 
sufficiently large solution to ip( z ) — 1 must be close to z m for some m, and since by construction there 
is only one such point z n inside the contour L(n), we conclude that for every n large enough there is a 
solution to ip(z) = q close to z n . □ 



Proof of Proposition 1: Let us assume that a < 1 and a < 1. Then we can take the cutoff function 
h(x) = in (2), and we can rewrite ip(z) as follows 



where we have denoted 



J (e zx - l)n(dx), 



(-00,0) 



{e zx - l)n(dx). 



(43) 

(44) 
(45) 



(0,*] 



First, let us study the asymptotic behavior of ip\{z). If a < then part (1) of Definition 2 implies 
that n((— 00, 0)) < 00, thus ipi(z) = 0(1) as z — > 00, Re(z) > 0. If < a < 1, then part (1) of Definition 
2 implies that for x G 1R~ 



n(dx) = ( Calx]- 1 ^ + J2^A\ X \ 



— l— &i 



dx + z/(dx), 



where u(dx) is a finite measure on K (note that z/(dx) does not have to be a positive measure). It is 
clear that 



{e zx - 1) u(dx) 



<2|i/((-oo,0))| 



for Re(,z) > 0. Using integration by parts we find that for Re(z) > 





a I (e zx -l)\x\- 1 - & dx = -T(l-a)z 



(46) 



Combining the above three equations and (44) we conclude that 

$ x {z) = -CT(1 - &)z & + o(z & ) 



(47) 
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as z — > oo, Re(z) > 0. 

Next, let us investigate the asymptotic behavior of ip2(z)- Let us assume that n > (where n is 
the constant in the Definition 2), the proof in the case n = is very similar. Since n > 0, Definition 2 
implies that the measure Il(dx) restricted to 1R + has a density it(x), which belongs to VC n ((0, k}). Let 
us assume first that it G C n ((0,k}), we will relax this assumption later. 

First let us consider the case a < 0, which is equivalent to n((0, oo)) < oo. Applying integration by 
parts n times to (44) we obtain 

k k 


n— 1 ^ ^ 

= JZ(-l) m [7r (m) (^-)e fc ^~ m_1 - 7r (m) (0+)^- m - 1 ] + (-l) n z~ n / 7r (n) (x)e^dx - / 7r(a;)da;. 

m=0 „ n 



Since 7r^ n \x) is continuous we conclude that 

k 

/^We-d, = .(«*) 
o 

as 2; — >■ 00, Re(^) > 0. At the same time, due to Definition 2 we have 7P m '(fc— ) = for all m < n — 2. 
Using the above two results and the fact that ^U + (x) = —tt(x) we obtain 

fo{z) = (-l) n Il +{n \k-)e kz z- n + o (e kz z' n ) + 0(1), (48) 

as z — > 00, Ke(z) > 0. Equation (48) shows that the exponential term in the right-hand side of (6) comes 
from the upper boundary of the support of the Levy measure and from the first non-zero derivative of 
n+(x) at k-. 

Next, let us assume that a 6 (0, 1). Then, according Definition 2, the density of the Levy measure 
can be expressed as follows 

m 

tt(x) = Cax^ a + C i a 3 x ~ a] + 9( x )> 

3=1 

where g G C ra ([0, k}). We can rewrite 4>2{z) as 

ip 2 (z) = CF(a,k,z) + '^2c j F(a j ,k,z) + / e zx g(x)dx - / g(x)dx, 

.7=1 { I 



(49) 



where we have defined 



F(a,M = a / (e --l),— d,. 



Let us obtain an asymptotic expansion of F(a, k, z) as z — > 00, z G Q\. Expanding exp(zx) in Taylor 
series centered at zero and integrating term by term we find that 

F(a, k, z) = k~ a [1 - iFnt-a, 1 - a; kz)] , (50) 
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where i-F\(a, 6; z) is the confluent hypergeometric function defined by (22). Applying asymptotic formula 
(2) on page 278 in [11] we conclude that 

F(a,k,z) = -r(l - a)e-™ z a + ak^— ]P 1 + ^ m (51) 

m=0 ^ ' 

+ 0(l) + 0(z a - l )+0(e z z~ N - 2 ), 

as 2 — >• oo, z G Qi. Formula (51) and our previous result (48) imply that 

fo( z ) = (-l)"n +(n) (A;-)e fe ^ n -r(l -a)e-™z a (52) 
+ o(e kz z- n ) +o{z a ) +0(1), 

as z — > oo, 2! G Qi. 

As a final step, let us relax the assumption ir G C n ([0, k]). Assume that there is a unique point 
x\ G (0, k) at which ix^ n \x) does not exist (the proof in the general case is exactly the same). According 
to Definition 2, tt G C n-2 (1R + ), thus 7r^(xi— ) = 7r( m )(:r 1 +) for m < n — 2. Applying integration by 
parts n times on each subinterval (0,o;i) and (xi,k) we would obtain an expression (52) plus an extra 
term of the form 

h(z) = (-I)"" 1 [vr^On-) - 7r (n - 1) (xi+)] e XlZ z- n . 

However, it is easy to see that h(z) = o(e kz z~ n ) +o(z a ) as z — > oo, z G Q%. This is true since in the domain 
V = {Re(z) < ln(ln(Im(z))), lm(z) > e} we have |exp(a;iz)| = exp(xiRe(2;)) = 0(ln|z|) = o(z a ), while 
in the domain Qi \ T> we have Re(z) — > oo when z — > oo, which implies exp(xiz)z~ n = o(exp(kz)z~ n ). 

Formula (43) and asymptotic expressions (52), (47) imply that ip(z) satisfies (6) with coefficients A, 
a, B and b as in Proposition 1, except that there would be an extra term 0(1) in the right-hand side 
of (6), which comes from (52). According to our assumption, the process X is not a compound Poisson 
process, thus the constant h defined in Proposition 1 is strictly positive, and 0(1) = o(z b ), therefore this 
extra term can be absorbed into o(z b ). This ends the proof in the case a < 1 and a < 1. 

In the case when one or both of a, a are greater than one the proof is identical, except that we will 
have to do one extra integration by parts for proving (47). The details are left to the reader. □ 



Proof of Theorem 3: Let us denote 



1 

k 



In 



+ (a + b) In 



2nir 



+ 



k 



arg 



+ \-(a + b) + 2n + 1 



7T 



Due to Definition 2, the Levy measure Il(da;) can only have a finite number of atoms. From Corollary 
2.5 in [23] we find that W^ q \x) can only have a finite number of points where it is not differentiable. 
Thus we can use (15) and the Bromwich integral formula to conclude that for any c > $(g) 

1 f zx 

dz. (53) 



x 



2711 



1py(z) - q 



For n > and m < we define the contour L = L(n, m) = L\ U L 2 U L 3 U L 4 , where 

Li = Li{n) = {Re(z) = c, — lm(z n ) — n/k < Im(z) < lm(z n ) + rr/k}, 

L 2 = L 2 (n,m) = {Im(z) = Im(z n ) + 7i/k, m < Re(z) < c}, 

L3 = Ls(n,m) = {Ke(z) = m, — lm(z n ) — n/k < lm(z) < Im(z„) + 7i/k}, 

L4 = L 4 (ra,m) = {Im(z) = — lm(z n ) —Tr/k, m < Re(z) < c}. 
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This contour is shown on figure 4b. We assume that L is oriented counter-clockwise. Using the residue 
theorem we deduce 



ip Y (z) 



-dz 



3 *(g)as 



-Coa; 



q 



+ 



+ 2^Re 



(54) 



where the summation is over all j > 1, such that — (j lie inside the contour L. 

First, assume that n is fixed and let us consider what happens as m — > —oo. According to the 
asymptotic relation (16), as Re(z) — > —oo the function ipy(z) increases exponentially (uniformly in every 
horizontal strip |Im(z)| < C). In particular, for m large enough we would have \iPy{ z ) — ^ | > 1 for all 
z G L^(n,m), which implies 



l[}y(z) 



-dz 



L 3 (n,m) 



Q 





e zx 




< I 




x \dz\ < j 




ipy(z) - q 



\dz\< J e xm x \dz\ = (2Im(z n ) + 2n/k) e r ' 



L 3 (n,m) 



L 3 (n,m) 



and for every x > the right-hand side converges to zero as m — > — oo. 

Our next goal is to let n — > +oo and to prove that the integrals over the two horizontal half-lines 
L 2 (n, — oo) and L^(n, — oo) in (54) disappear. In order to achieve this we'll need to obtain good upper 
bounds on ip(z) on these horizontal half-lines. Let us consider first the contour L 2 (n, — oo). We will 
prove that there exists a constant C such that |-?/>y(z)| > C|hm(z n )| for all z G L,2(n, — oo). 

Assume that e > is a small number and define a domain 



V 6 = {z G C : |arg(z)| > vr/2 + e}, 

see figure 4b. Let L 5 = L 5 (n) = L 2 {n ) —oo) fl V e and L 6 = L e (n) = L 2 (n, — oo) \ V e . Following the same 
steps as in the proof of Theorem 2 (see estimate (42)) we find that there exists a constant C\ (which 
does not depend on n or e) such that for all n large enough we have for all z G L^(n) 

\ifr(z)\ > cos(Cie) (\Ae kz z~ a \ + \Bz b \) > cos{C 1 e)\B\\z\ b (55) 
> cos(Cie)|fi|Im(z n ) b > cos(Cie)| J B|Im(z n ) 

where in the last estimate we have used the fact that b > 1 (see (17)). 

Next, it can be easily seen from the figure 4b that for all z in the domain T> e we have 

Re(^) < — \z\ sin(e), (56) 

therefore \z\ b = o(exp(—kz)z~ a ) when z — > oo, z G T> e . This fact and the asymptotic formula (16) show 
that there exists a constant C 2 such that for all z G V t large enough we have |-0y(,z)| > C 2 \ exp(— kz)z~ a \. 
Therefore, for all n large enough we have 

\ij Y (z)\>C 2 e klKe{z)l \z\- a , zeL 6 (n). 

Using the above estimate and (56) we find 

\*Py{z)\> C 2 sm{e) a e klRe(z)l \Re{z)\- a , z G L 6 (n). (57) 

Next, as n increases to +oo, the real part of any z G L 6 (n) decreases to — oo (see figure 4b), thus for all 
n large enough we have exp(/c|Re(;?)|/2) > |Re(z)| a for all z G Lq(u). At the same time, from the figure 



25 



4b we see that for all z e L§{n) it is true that |Re(z)| > tan(e)|Im(z)| = tan(e)(Im(^ n ) +ir/k). Using 
this fact and (57) we find that there exists a constant C3 = C^{e) such that for all n large enough 

\ifr(z)\ > C 2 sin(e) a etl Re ( 2 )l|Re(z)r a e5l Re ( 2 )l > C 2 sm(e) a e^ tffl(f),ffi(z " ) > C 3 Im(z n ), z e L 6 (n). (58) 

Combining (55) and (58) we conclude that there exists a constant C > 0, such that for all n large enough 
we have 

\i/) Y (z)\ > CIm(z n ), z G L 2 (n, —00). 

A similar estimate for L^(n, —00) can be obtained in the same way. 
Thus setting z = z(u) := u + i(Im(;z n ) + n/k) we obtain 



1,2(11,— 00) 

c 



ip Y (z) -q 



dz 



z(u)x 



du 



< 



7 z(u)x I 



\^ Y {z{u)) - q 



-du < 



-du 



C\Im{z n )\-q ■ C\lm(z n )\-q 



and the right hand side converges to zero asn-> +00. Similarly, the integral over L 4 vanishes. Thus as 
n — > +00 formula (54) becomes 



-dz 



c+i]E 



n>l 



and the left-hand side is equal to W^ q '(x) due to Bromwich integral formula (53). 

Next, from Proposition 1 we know that the asymptotic formula for ip Y ( z ) can be obtained by differ- 
entiating (16). Therefore, using the asymptotic expression (7) for ( n we find that 



M-Cn)\=k\B\ 



2mr 



+ o (n b ) 



Similarly, from (7) we find that there exists a constant c such that 



\e <nX \ ~ cn fc (a+,,) , n -> +00, 



thus the terms of the series in the right-hand side of (18) decrease as n - b - x ( a + b )/ k , According to (17) 
we always have b > 1, which implies that the series in the right-hand side of (18) converges on M + and 
uniformly on [e, 00) for each e > 0. □ 

Proof of Proposition 2: First, let us assume that a < 1 and a < 1. We start with the Levy-Khintchine 
formula (2) with the cutoff function h(x) = 0, which gives us 



il>(z) = -a 2 z 2 + iiz + Ca I (e zx - 1) e^'^dx + Ca \ (e zx - 1) e^x'^dx. (50) 
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The first integral in (59) can be evaluated as follows: 

oo o 
Ca J (e zx -l)e^ x \x\^- & dx = Ca J (e {z+ ? )x - l) {x^^dx - Ca J [eP x - l) \x\~ l ~ & dx 



oo 

a. 



= -CT(1 -a)0 + zf + CT(1 -&)P 

where we have used (46) in the final step. Similarly, the second integral in (59) can be evaluated with 
the help of (21) and (50): 

k k k 

Ca J (e zx - 1) e-^x-^dx = Ca J ( e ^~ p)x - l) x'^dx -Ca J (e~ Px - l) x'^dx 



o 



= CF(a,k,z-P)-CF(a,k,-P) 

= -Ck^F^-a, 1 - a, -k(/3 - z)) + Ck~ a x F x {-a, 1 - a, -hp) 
= Ca(f3 - z) a <y(-a } k(P - z)) - Ca(3 Q -f(-a, k/3). 

This ends the proof in the case a < 1 and a < 1. When a > 1 or a > 1 the proof would be very 
similar, the only difference is that we would use the cutoff function h(x) = 1 in (2) and perform an extra 
integration by parts in (46). We leave all the details to the reader. □ 
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